Resonant Einstein-de Haas effect in a rubidium condensate 
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We numerically investigate a condensate of *^Rb atoms in an _F = 1 hyperfine state confined in 
an optical dipole trap. Assuming the magnetic moments of all atoms are initially aligned along 
the magnetic field we observe, after the field's direction is reversed, a transfer of atoms to other 
Zeeman states. Such transfer is allowed by the dipolar interaction which couples the spin and the 
orbital degrees of freedom. Therefore, the atoms in mp = 0,-1 states acquire an orbital angular 
momentum and start to circulate around the center of the trap. This is a realization of the Einstein- 
de Haas effect in systems of cold gases. We find resonances which amplify this phenomenon making 
it observable even in very weak dipolar systems. The resonances occur when the Zeeman energy on 
transfer of atoms to mp = state is fully converted to the rotational kinetic energy. 
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Magnetic effects in ultracold quantum gases have been 
subject to intense theoretical and experimental studies 
during recent years. So far most of these investiga- 
tions have concentrated on short-range interactions as 
the dominant spin exchange process in spinor conden- 
sates of 23Na [l| and ^^Rb [1, i, i]. These interac- 
tions result in rich multi-component physics as demon- 
strated by the observation of phenomena like magnetic 
phases [II H, 0] ) coherent spin dynamics 0, Hi S , do- 
main formation [3] , and a magnetically tuned spin mixing 
resonance Q. Including magnetic dipole-dipole interac- 
tions would further enhance the richness of these systems 
and in particular their anisotropic nature is expected to 
add completely new aspects. Already for relatively weak 
dipolar interactions phenomena like the Einstein-de Haas 
effect [§] , spontaneous magnetization [10,, JT,] , squeezing 
and entanglement [ill have been predicted. Most of these 
studies concentrate on the recently achieved case of a 
chromium Bose-Einstein condensate [l2| . as it is com- 
monly believed that these effects are practically unob- 
servable in the widely available alkali condensates due 
to the smallness of dipolar interactions in these systems. 
However, as pointed out in [111, for ^"^Rb in the F = 1 
state the size of the dipolar interactions as compared 
to the spin-mixing part of the short-range interactions 
reaches 10%, such that dipolar effects might be observ- 
able in this system. 

In this Letter, we show that under the right conditions 
the dipolar interactions can even dominate the dynamics 
of a *^Rb spinor condensate, making it a promising can- 
didate for the observation of the Einstein-de Haas effect 
[T^ as an important milestone in dipolar quantum gas 
physics. In particular our calculations demonstrate the 
existence of resonances in a spinor rubidium condensate 
that amplify the effect of dipolar interactions. These res- 
onances can be tuned by the magnetic field or by the trap 
geometry and occur when the Zeeman energy fits the ro- 
tational kinetic energy per particle. The resonances we 
find explore a new regime in comparison with that con- 
sidered in Ref. [H for a ^^Cr condensate, where the dipo- 
lar energy (not a kinetic one) is related to the Zeeman 



energy. 

In the second quantization notation, the Hamiltonian 
of the system we investigate is given by 
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where repeated indices (each of them going through the 
values +1, 0, and —1) are to be summed over. The field 
operator (r) annihilates an atom in the hyperfine state 
\F = l,z > at point r. The first term in ^ is the 
single-particle Hamiltonian (7?o) that consists of the ki- 
netic energy part (with m being the mass of an atom) 
and the trapping potential Uext (r) . The second term de- 
scribes the interaction with the magnetic field B with 7 
being the gyromagnetic coefiicient which relates the ef- 
fective magnetic moment with the hyperfine angular mo- 
mentum (/X = 7F) . The terms with coefficients cq and C2 
describe the spin-independent and spin-dependent parts 
of the contact interactions, respectively - cg and C2 can be 
expressed with the help of the scattering lengths oq and 
a2 which determine the collision of atoms in a channel of 
total spin and 2. One has cq — 47rfi.^(ao + 202) /3m and 
C2 = 47r/i^(a2 — ao)/3m [l^, where ao = 5.387 nm and 
02 = 5.313nm 15. F are the spin-1 matrices. Finally, 
the last term describes the magnetic dipolar interactions. 
The interaction energy of two magnetic dipole moments 
IJ,i and ^2 positioned at r and r' equals 
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and since /x = 7F one has Vf,j.i{Y - r') = 7^FjjFfe;/|r - 
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The equation of motion reads 

d ft' \ (i^ \ 

'^di ^" j=i^- + ^d) 1^ ^0 j , (3) 

where the operator Tic originates from the single-particle 
and two-particle but related to the contact interactions 
parts of the Hamiltonian ^ whereas TLd corresponds to 
the dipole-dipole interactions. The diagonal part of Tic 
is given by Hcii = Hq + (cq + C2) V'lV'i + (co + C2) i/'(5V'o + 

(co-C2)'0li^/'-l, UcOO = i7o + (co + C2)i/'li/'i-f Co'(/'o?/'0 + 

(co -I- C2)^/'li^/'-i, Uc-i-i = Ho + (cq - C2) + (co + 

C2) V'oV'o + (co + C2) The off-diagonal terms that 

describe the collisions not preserving the projection of 
spin of each atom (although the total spin projection is 

conserved) equal Hcia = C2'4>''_i^o, Wco-i — C2'4'o'4'i- 
Moreover, Hci-i = 0. On the other hand, for the Hd 
term one has Hdij — f d^r''^lJ^^{r')V^j ^i,^k{^')- This term 
is responsible for the change of total spin projection of 
colliding atoms. It turns out that when two atoms inter- 
act the total spin projection {AMp) can change at most 
by 2. In particular, the diagonal elements of Hd lead to 
the processes with AMp = ±1. An example could be 
the collision of two atoms in m^? = 1 Zeeman state after 
which one of the atoms goes to the component mp — 0. 
In addition to such processes, the off-diagonal terms of 
Tid introduce the interaction that change the spin projec- 
tion AMp by ±2. It happens when both atoms initially 
in the same state go simultaneously to the nearest (in a 
sense of magnetic number mp) state or in the case when 
atoms in different but neighboring components transfer 
to the states shifted in number mp by -1-1 or —1. There 
is no way for the atom being transfered directly from 
mp = 1 to the mp = —1 state, therefore, the populating 
of mi? — —1 component is a second order process. 

Hence, the dipolar interaction does not conserve the 
projection of total spin of two interacting atoms. Nei- 
ther the projection of total orbital angular momentum is 
preserved (see Q). However, the dipolar interaction cou- 
ples the spin and the orbital motion of atoms as revealed 
by the last relation in ^ 

[Vd, Fi, + ^ 0, [Vd,Li, + L2,]^0, 

[Vd, Li, + L2, + Fi, + F2,] ^ . (4) 

Therefore, going to mp = 0,-1 states atoms acquire the 
orbital angular momentum and start to circulate around 
the center of the trap. This is the realization of the fa- 
mous Einstein-de Haas effect in cold gases. 

To solve the Eq. ^ we neglect the quantum fluctu- 
ations and replace the field operator "01 (r) by an order 
parameter ipi (r) for each component and apply the split- 
operator method. All integrals appearing in Tidij are 
the convolutions and we use the Fourier transform tech- 
nique to calculate them. To find analytical formulas for 
the Fourier transforms of the components of the convo- 
lutions that do not change during the evolution we apply 
the regularization procedure described in Ref. p^ . 



The gyromagnetic coefBcient for Rb atoms in an _F = 
1 hyperfine state is positive and equals 7 = i/Lts/fi. We 
prepare an initial state of the condensate as the one with 
all magnetic moments aligned along the magnetic field, 
i.e., all atoms are in mp — 1 component. To this end, 
we run the mean-field version of Eq. ([3|) in imaginary 
time while the magnetic field is turned on (and equal 
to B = 0.73 mG for a spherically symmetric trap with 
the frequency tu = 2tt x 100 Hz). Then we reverse the 
direction of the magnetic field and look for the transfer 
of atoms to other Zeeman states. 
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FIG. 1: (color online). Transfer of atoms to mp ~ 0,-1 
Zeeman states as a function of time. Initially, N = 2 x 10^ 
atoms were prepared in = 1 component in a spherically 
symmetric trap with the frequency uj = 2n x 100 Hz. The 
residual magnetic field equals —0.029 mG (black lines - on- 
resonance case) and B = —0.015 mG, —0.036 mG (blue and 
red lines, respectively - ofF-resonance case). The lower panel 
shows: the number of atoms in mp — 1 state and the time 
dependence of the orbital angular momentum per atom. 

Our starting condition (all atoms in mp = 1 state) 
suppresses the short-range spin dynamics and initially 
the mp = \ state is depleted only due to the dipolar in- 
teraction. Usually, we observe a small number of atoms 
going from the Zeeman state mp — \ to the mj? = 0, — 1 
states. However, on resonance (see Fig. [1} the transfer 
to the other states can be of the order of the initial popu- 
lation oi mp = 1 component. This transfer is as large as 
in the case of chromium condensate @ despite the fact 
that the dipolar energy (/x^n, where n is the atomic den- 
sity) is approximately 100 times smaller. The only differ- 
ence is that the time scale corresponding to the maximal 
transfer is about 100ms, i.e., 100 times longer than for 
chromium. This can be understood as follows. For *^Rb 
the dipolar energy is the smallest energy of the problem 
and is only a perturbation as compared to the kinetic 
rotational energy and the Zeeman energy. An efficient 
population transfer from mp = 1 to mp — Q due to the 
dipolar interactions is only possible, when the total en- 
ergies (the sum of mean-field, trap, kinetic, and Zeeman 
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energies) in these states are approximately equal. Fig. 
[5^ clearly shows that the total energy tends to equalize 
on resonance, which is not true in the off-resonant case. 
Numerics also shows (Fig. that the above condition 
can be fulfilled only then when the Zeeman energy fits 
the rotational kinetic energy: ^Bres = f-rot- ft means 
that the resonant magnetic field is inversely proportional 
to the magnetic moment of an atom 



(5) 



Surprisingly, the smaller atomic magnetic moment the 
larger value of the resonant magnetic field. For chromium 
condensate, however, the dipolar energy is larger than 
the kinetic energy. Therefore, in this case the resonance 
condition should be derived by relating the Zeeman and 
the dipolar (not the kinetic one) energies, i.e., ^B^es — 
fi^n [9]. This results in a condition Bres oc M which differs 
qualitatively from ([5|). 

The maximal transfer is reached at time which is of 
the order of characteristic time scale determined by the 
dipolar interactions (fi/ fj?n). Since the magnetic moment 
of ^^Rb atom is 12 times smaller than that of ^^Cr we 
have to wait hundreds of milliseconds (not a fraction of 
millisecond as in Ref. Q) to see the action of resonance. 

Fig. [5] illustrates the ideas just discussed. In the up- 
per frame the total energies oi mp — 1,0 components 
are plotted as a function of time both in on- and off- 
resonance cases showing that the resonances we find is 
a dynamical phenomenon. The transfer gets maximal 
when the energies approach each other (perhaps crossing 
both curves would require the dynamical tuning of the 
resonance by changing the magnetic field). Contrary, al- 
most no transfer of atoms occurs when the energy curves 
keep away. Simultaneously, the lower frame proves that 
on resonance the Zeeman energy (in fact, together with 
the kinetic energy) is transfered to the rotational kinetic 
energy of atoms in mp = component. 

Huge transfer of atoms to mp — 0,-1 states is the 
realization of the Einstein-de Haas effect in cold gases. 
Fig. [5^ proves that the vortices are generated in m^^ = 
0,-1 components (though already the lower-right frame 
in Fig. [l]suggests that atoms in states nip = 0, —1 rotate 
around the quantization axis). For mp = component 
the phase of the order parameter winds up by 27r which 
means that a singly quantized vortex is formed in this 
state. At the same time in the mp = —1 component a 
doubly quantized vortex is formed (the phase winds up by 
47r) as a result of total angular momentum conservation. 
Fig. shows the typical density patterns in the mp ~ 
0,-1 states. The density is fragmented and the number 
of rings in the mp = 0,-1 components results from the 
symmetry of dipolar interaction. For example, since the 
dipolar interaction transforms as a spherical tensor of 
rank-2 in both spatial and spin spaces, it spatially acts as 
the spherical harmonic 1^21 when atoms go from mp = 1 
to mp = state. And, since the initial state is spherically 
symmetric (oc loo) it induces the spatial behavior oc 
in mi? = component. Therefore, the density inmp = 
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FIG. 2: (color online). Total energy for mp = 1,0 com- 
ponents as a function of time (upper frame). All param- 
eters are the same as in Fig. [1] The resonance happens 
for B = —0.029 mG (black lines) whereas the off-resonance 
cases are represented hy B = —0.015 mG (blue lines) and 
B — — 0.036 mG (red lines). The lower frame shows AE = 
{E^l^/N+i + fiB - erot)/€rot, wherc erot = Erot/No. 




FIG. 3: (a) Phase in the 'xy' plane of ttif = (left frame) 
and jTiF ~ —1 (right frame) spin components, (b) Density 
in the 'xz' plane (the 'z' axis goes vertically). AH parameters 
are as in Fig. [T] (on-resonance case) and the density cuts are 
taken at 140 ms. 



state shows two rings. Similar fragmentation was already 
predicted in the case of ^^Cr condensate in Ref. Q. 

Fig. m (upper frame) shows the position and the width 
of the resonance displayed in Fig [TJ Similar behavior 
is observed when the value of the reversed magnetic 
field is kept constant and the trap geometry is changed 
(lower frame). Here, the maximal transfer of atoms is 
obtained in a cigar trap with the aspect ratio tUp/tUz = 4 
with almost 50% efficiency at B = -0.073 mG. The in- 
set shows the resonance at experimentally easier to con- 
trol value of magnetic field B = 0.3 mG but still de- 
tectable number of atoms in m^? = state. To un- 
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derstand quantitatively the resonance we start from the 
condition discussed ear her: /ii?res = Erot/No, where 
Erot is the rotational energy of = component 
which is assumed to be a singly quantized vortex, given 
within the Thomas- Fermi approximation by ipoiPi 4>i ^) = 
[(A - mtj2(p2 _^ - j {2mp^)) I c^fl"^ x e"^. Here, 

the chemical potential A is obtained by the requirement 
that the number of atoms in state equals -/Vq. One 
can tune to the resonance in two ways: a) by adjusting 
the magnetic field B\ b) by changing the trap geometry 
what influences the rotational energy entering resonance 
condition and keeping the magnetic field constant. The 
resulting curve for the spherically symmetric trap is plot- 
ted in Fig. [5l The numerical results are marked by balls 
and correspond to the initial number of atoms in tof = 1 
component equal to 2 x 10^, 8 x 10^, and 2 x 10^ showing 
a good agreement. For other systems, e.g. ^^Cr, the con- 
dition iiBres — Erot/Nmp suggcsts that the value of the 
resonant magnetic field is even k, 10 times smaller since 
fiCr/fJ-Rb = 12 and Erot /N -2 for chromium looks similar 
as Erot 1^0 in the rubidium case. 
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FIG. 4: Maximal transfer of atoms to mp = 0,-1 states 
as a function of the residual magnetic field (upper frame) 
and the trap geometry (lower frame). For lower frame lu^ — 
2tv X 100 Hz, B = -0.073 mG, and A^+i = 8 x 10^ Inset 
shows the resonance at _B = —0.3 mG for ujz = 2n x 20 Hz 
and N+i = 10^ 



In conclusion, we have shown the existence of dipo- 
lar resonances in rubidium spinor condensates. The 
resonances occur when the Zeeman energy of atoms in 
mp — 1 component, while transferring to mp = state, 
is fully converted to the rotational kinetic energy. This 
is so far an unexplored regime. Symmetries of the dipo- 
lar interaction force the atoms in mp = 0,-1 states to 
circulate around the quantization axis and form singly 
and doubly quantized vortices, respectively. Therefore, 
dipolar resonances is a route to the observation of the 
Einstcin-dc Haas effect (as well as other phenomena re- 
0.06 
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FIG. 5: Comparison between numerics (balls) and the 
Thomas-Fermi approximation. Solid line indicates the value 
of the magnetic field at resonance (for spherically symmetric 
trap) as a function of number of atoms in mp =0 state. 



lated to the dipolar interaction) in weak dipolar systems. 
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